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We study the propagation of charm quarks, produced from 
the initial fusion of partons, in an equilibrating quark-gluon 
plasma which may be formed in the wake of relativistic col¬ 
lisions of gold nuclei. Initial conditions are taken from a self 
screened parton cascade model and the chemical equilibration 
is assumed to proceed via gluon multiplication and quark pro¬ 
duction. The energy loss suffered by the charm quarks is ob¬ 
tained by evaluating the drag force generated by the scatter¬ 
ing with quarks and gluons in the medium. We find that the 
charm quarks may loose only about 10% of their initial energy 
in conditions likely to be attained at the Relativistic Heavy 
Ion Collider, while they may loose up to 40% of their energy 
while propagating in the plasma created at Large Hadron Col¬ 
lider. We discuss the implications for signals of quark gluon 
plasma. 

PACS numbers: 12.38.Mh, 12.38.-t, 14.65.Dw, 25.75.-q 


I. INTRODUCTION 

Relativistic heavy ion collisions are being studied with 
the intention of investigating the properties of ultradense, 
strongly interacting matter- quark-gluon plasma (QGP) 
[|ij. The last several years have witnessed extensive the¬ 
oretical efforts towards modeling these collisions. These 
studies suggest that at collider energies, one may visual¬ 
ize the nuclei as two clouds of valence and sea partons 
which pass through each other and interact [p|. The re¬ 
sulting partonic interactions may then produce a dense 
plasma of quarks and gluons. This plasma will expand 
and become cooler and more dilute. If quantum chro¬ 
modynamics admits a first-order deconfinement or chiral 
phase transition, this plasma will pass through a mixed 
phase of quarks, gluons, and hadrons, before the hadrons 
loose thermal contact and stream freely towards the de¬ 
tectors. 

One would like to understand several aspects of this 
evolution, viz., how does the initial partonic system 
evolve and how quickly does it attain kinetic equilibrium? 
How quickly, if at all, does it attain chemical equilibrium? 
And finally, how can we uncover the history of this evo¬ 
lution by studying the spectra of the produced particles, 
many of which may decouple from the interacting sys¬ 
tem only towards the end? These and related questions 
have been actively debated in recent times. Thus, it is 
believed by now that the large initial parton density may 
force many collisions among the partons in a very short 
time and lead to a kinetic equilibrium Q. The question 
of chemical equilibration is more involved, as it depends 


on the time available to the system. If the time available 
is too short (3-5 fm/c), as at the energies (-^s < 100 
GeV/nucleon) accessible at the Relativistic Heavy Ion 
Collider (RHIC), the QGP will end its journey far away 
from chemical equilibrium. At the energies (y/s < 3 
TeV/nucleon) that will be achieved at the CERN Large 
Hadron Collider (LHC) this time could be large (more 
than 10 fm/c), driving the system very close to chemical 
equilibration (J |fj, if there is only a longitudinal expan¬ 
sion. However, this time is also large enough to enable a 
rarefaction wave from the surface of the plasma to prop¬ 
agate to the center (r s ~ Rt/c s ; Rt is the transverse 
nuclear dimension and c s is the speed of sound), thus 
introducing large transverse velocity (gradients) in the 
system. The large transverse velocities may impede the 
chemical equilibration by introducing a faster cooling. 
The large velocity gradients may drive the system away 
from chemical equilibration by introducing an additional 
source of depletion of the number of partons Q in a given 
fluid element. 

Can we identify a probe for these rich details of the 
evolution? Dileptons and single photons have long been 
considered as useful probes of the various stages of the 
plasma as once produced they hardly ever interact and 
thus carry the details of the circumstances of their pro¬ 
duction. They are however produced at every stage of 
the collision and in an expanding system their number is 
obtained by an integration over the four-volume of the 
interaction zone. This moderates our estimate of the ef¬ 
ficacy of these signals in several competing ways, as the 
space-time occupied by the hot plasma is small. 

Consider the case of single photons. At very early 
times the temperature is rather high and we have a co¬ 
pious production of photons having a large transverse 
momentum. This should give us a reliable information 
about the initial conditions of the plasma. However, the 
transverse flow of the system is very moderate at early 
times. By the time the flow and other aspects of the 
QGP develop, the temperature would have dropped con¬ 
siderably and we have a large production of photons hav¬ 
ing smaller transverse momenta. However, the photons 
having a small transverse momentum will also be pro¬ 
duced, and in much larger numbers, from the hadronic 
reactions and decays of vector mesons in the hadronic 
matter. Thus disentangling the development of the later 
stages of the evolution in the QGP would entail a very de¬ 
tailed understanding of the contributions of the hadronic 
matter. This is not trivial. 

Next consider dileptons. Again, large mass dileptons, 
having their origin in early times, will be only marginally 
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affected by these developments, e.g, the flow. The low- 
mass dileptons produced at later times will be affected 
strongly by the flow. However they will have, in addition 
to the contributions from the plasma, a large contribu¬ 
tion from hadronic reactions. We shall be required to un¬ 
derstand them in detail before we can confidently embark 
upon the task of deciphering the, more involved, develop¬ 
ment of the QGP, as it gets cool (see, e.g.,Ref, |J). At the 
same time, there have been suggestions that correlated 
charm or bottom decay may give a large background to 
the dileptons from the QGP. We shall come back to this 
aspect later. 

Why should we be so concerned about these late de¬ 
velopments in the plasma? Firstly, it is not difficult to 
imagine that the detailed constitution of the hadronic 
matter should somehow reflect the state of chemical equi¬ 
libration of the QGP at the time of the phase transition. 
Secondly, the transverse velocities, which may develop 
during the QGP phase will remain unaltered during the 
mixed phase, if such a phase is attained. The develop¬ 
ment of the transverse velocities should also be affected 
by the speed of sound, which may change rapidly as we 
approach the transition temperature. In the absence of 
transverse velocities, the mixed phase will last for a very 
long period of time |il)| . And lastly, it is not at all clear 
that a non-equilibrated plasma will be adiabatically con¬ 
verted to hadronic matter, a scheme which is often in¬ 
voked in hydrodynamic descriptions of the system. 

In the present work we report on an investigation of 
the energy loss suffered by charm quarks in a chemically 
equilibrating plasma. This study has several interesting 
features. Firstly, heavy quarks are mostly produced from 
the initial fusion of partons in the colliding nucleons. If 
the initial temperature is high, at least charm quarks may 
also be produced at very early times (mostly from gg —> 
cc, but also from qq —» cc). There is no production of 
charm quarks at later times and none in the hadronic 
matter. This makes them a good candidate for a probe 
of QGP. Again, as their number is not large, one can 
confidently ignore the reverse process (cc —> ggorqq), 
where q denotes one of the lighter quarks. (Of-course the 
associated production/suppression of J/’F and T is a sub¬ 
ject by itself.) The number of initially produced charm 
quarks (from fusion of primary partons) can be evaluated 
with a great degree of confidence using the perturbative 
QCD |nj. The thermally produced heavy quarks can 
also be estimated in terms of the temperatures and the 
fugacities during the early stage of evolution, when they 
are produced. Thus, the total number of c or b quarks 
gets frozen very early in the history of the collision. This 
can be of great help, as we are then left with the task 
of determining their px distribution, whose details may 
reflect the developments in the plasma, but which is nor¬ 
malized to a given number of c (or b) quarks. 

Now consider the central slice (y = 0) of a collision 
which leads to production of a charm quark pair from ini¬ 
tial fusion and QGP. The charm quarks will be produced 
on the order of a time scale 1/2 to c ~ 0.07 fm/c, while the 


bottom quarks will similarly be produced over a time of 
~ 0.02 fm/c. Thus, immediately upon production, these 
quarks will find themselves in a deconhned matter which 
is rapidly evolving, first towards kinetic equilibration and 
then towards chemical equilibration. The heavy quarks, 
as well as the light quarks and gluons will undergo scat¬ 
tering, with an interesting difference. The heavy quarks 
are not very likely to have a frequent and drastic change 
in their direction upon bombardment by massless quarks 
and gluons. These collisions will amount to a drag force. 
Svetitsky iQ has argued that this drag force is rather 
large and it will bring down the velocity of charm quarks 
down to thermal velocity very quickly. 

In the present work we show that when we account 
for the non-equilibrium nature of the plasma, this drag 
force reduces considerably. In particular, we find that the 
charm quarks produced at the RHIC may loose only up to 
10% of their initial momentum though at the LHC they 
may loose up to 40% of their initial momentum during 
their passage through the QGP phase of the matter. 

The momenta of these quarks are likely to be reflected 
in the corresponding quantities of D mesons as the c 
quarks should pick up a light quark, which are in great 
abundance and hadronize. (One may add an interesting 
scenario, where the D mesons travel through a string of 
QGP droplets and loose further energy |l3|].) Finally, 
they would scatter with hadronic matter before decou¬ 
pling. Considering that cjd-k "C cr W7r , e.g., it is most likely 
that the heavy mesons would decouple quickly from the 
hadronic phase. 

There has been an attempt | |i~i}[K| to estimate the 
change in the distribution of D mesons (open charm) by 
introducing a rate of energy loss dE/dx « —2 GeV/fm for 
the charm and bottom quarks. We recall that this value 
was originally obtained [fl6| for very energetic quarks. 
For this reason, this estimate also derives a considerable 
contribution from the radiation of gluons. Charm quarks 
are very massive, and thus the radiation mechanism is 
absent 0 unless their energy is exceptionally high. We 
shall be concerned with charm quarks having a px of a 
few GeV/c, in our study. For these energies, the drag 
force will come mostly from scatterings with light quarks 
and gluons |l2| in a formulation which is equivalent to 
the treatment of Braaten and Thoma [Q for this pur¬ 
pose. We show in the present work that this energy loss 
can be very low for a chemically non-equilibrated plasma. 

In the next Section, we briefly recall the initial condi¬ 
tions, and hydrodynamic and chemical evolution of the 
plasma in a (1+1) dimensional expansion. In Section III 
we give the formulation for drag force operating on the 
charm quarks in such an equilibrating plasma. We dis¬ 
cuss our results in section IV. Finally, in section V we 
give a brief summary. 
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II. INITIAL CONDITIONS, HYDRODYNAMIC 
EXPANSION, AND CHEMICAL 
EQUILIBRATION 

It is quite clear that the fate of the heavy quarks will 
depend on the initial conditions and the history of evolu¬ 
tion of the plasma. Fortunately, by now, a considerable 
progress has been achieved in our understanding of the 
parton cascades which develop in the wake of the colli¬ 
sions. Thus, in the recently formulated self-screened par- 
ton cascade |il| model early hard scatterings produce a 
medium which screens the longer ranged color fields as¬ 
sociated with softer interactions. When two heavy nuclei 
collide at sufficiently high energy, the screening occurs 
on a length scale where perturbative QCD still applies. 
This approach yields predictions for the initial conditions 
of the forming QGP without the need for any ad-hoc mo¬ 
mentum and virtuality cut-off parameters |2|]. These cal¬ 
culations also show that the QGP likely to be formed in 
such collisions could be very hot and initially far from 
chemical equilibrium. 

We assume that kinetic equilibrium has been achieved 
when the momenta of partons become locally isotropic. 
At the collider energies it has been estimated that, r* ss 
0.2 — 0.3 fm/c ||. Beyond this point, further expansion 
is described by hydrodynamic equations and the chemi¬ 
cal equilibration is governed by a set of master equations 
which are driven by the two-body reactions ( gg «-> qq) 
and gluon multiplication and its inverse process, gluon 
fusion (gg «-> ggg). The other (elastic) scatterings help 
maintain thermal equilibrium. The hot matter continues 
to expand and cools due to expansion and chemical equi¬ 
libration. We shall somewhat arbitrarily terminate the 
evolution once the energy density reaches some critical 
value (here taken as e/ = 1.45 GeV/frn 3 [pof ). 

The expansion of the system is described by the equa¬ 
tion for conservation of energy and momentum of an ideal 
fluid: 

d^ v = 0 , T» v = (e + P)u»u v + Pg^ , (1) 

where e is the energy density and P is the pressure mea¬ 
sured in the frame comoving with the fluid. The four- 
velocity vector of the fluid satisfies the constraint 
u 2 = —1. For a partially equilibrating plasma the dis¬ 
tribution functions for gluons and quarks can be scaled 
through equilibrium distributions as 


g i (q,T,X i ) = X i g^( q ,T) , (2) 

where g° q (q,T) = (e^ u ' q =p l) -1 is the BE (FD) dis¬ 
tributions for gluons (quarks), and A» is the fugacity 
for parton species i, which describes the deviation from 
chemical equilibrium. This fugacity factor takes into ac¬ 
count undersaturation of parton phase space density, i.e., 
0 < A* < 1. The equation of state for a partially equi¬ 
librated plasma of massless particles can be written as 

I 


e — 3 P — \a,2^g + (A q + Xq)] T 4 , 

where a 2 = 87t 2 /15, b 2 = 77r 2 TV//40, Nf ss 2.5 is the 
number of dynamical quark flavors. Now, the density of 
an equilibrating partonic system can be written as 

Tig — XgTlg, Tlq — Xqflq , (4) 

where hk is the equilibrium density for the parton species 
k: 

1 fi 

hg = — C(3)T 3 = ai T 3 , (5) 

TT 


n q = 2~2 C(3)A/T 3 = bi T 3 . (6) 

We further assume that \ q = Xq. The equation of state 
(||) implies the speed of sound c s = 1/ V3. 

The master equations [Q for the dominant chemical 
reactions gg <-> ggg and gg <-> qq are 


dg,(Tlg1l^') — Tlg(R2 —>3 P'3 —>2) ( TLgRg — TlqRq — >g) , 

dqfaqll 1 *) = dg(nqU^) = TLgRg —>g TlqRq —(7) 


in an obvious notation. 

If we assume the system to undergo a purely longitu¬ 
dinal boost invariant expansion, (jl) reduces to the well 
known relation 


de 

dr 


e + P 


= 0, 


where r is the proper time. This equation implies 
r 4 / 3 = const. 


e t 


and the chemical master equations reduce to 


1 dX g 
X g dr 
1 dX n 


3 dT 1 


T dr 


+ - = R 3 ( 1 - a g) - 2R,2 ( 1 - 


XnX, 


q''q 


% 


3 dT 1 


a i 


T-+ - = 7?2 7“ I T“ "T 2 " 


X q dr 


T dr 


X n X n 


bl V Xq Xq 


( 8 ) 

(9) 


( 10 ) 


which are then solved numerically for the fugacities with 
the initial conditions obtained from self screened parton 
cascade model We quote the rate constants i ?2 and 
I ?3 appearing in Eq. (0) for the sake of completeness [Q; 

R 2 « 0.2dN f a 2 s X g Thi(1.65/a s Xg), 

R 3 = 1.2a 2 T(2A 5 — A 2 ) 1 / 2 , (11) 

where the color Debye screening and the Landau - Pomer- 
anchuk - Migdal effect suppressing the induced gluon ra¬ 
diation have been taken into account, explicitly. 
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III. DIFFUSION OF HEAVY QUARKS IN 
QUARK GLUON PLASMA 


The early results of parton cascade model j|] do in¬ 
dicate that the charm quarks which are produced from 
initial fusion in heavy ion collisions may have a transverse 
momentum distribution given by a power law, initially. 
With the passage of time, these distributions start evolv¬ 
ing into an exponential shape, due to scattering with 
other partons, which is characteristic of a thermalized 
system of particles. Will these heavy quarks stay in ther¬ 
mal equilibrium, till the end of the QGP phase? For this 
we have to continue the cascading till the temperatures 
have dropped to T c or the energy density has become too 
low. The parton cascade model also includes all the ef¬ 
fects like scatterings and radiation, if any, on the heavy 
quark. We shall report on such an effort in a future pub¬ 
lication p^| . 

In the present work, we adopt the procedure developed 
by Svetitsky jT^] , where one visualizes the effect of par- 
tonic collisions as leading to a drag force. We first gen¬ 
eralize the treatment of Svetitsky to a non-equilibrated 
plasma and then evaluate the time variation of the so- 
called drag and the diffusion coefficients, to see whether 
heavy quarks will actually stop and diffuse at RHIC and 
LHC energies. 

We write the Boltzmann equation for the density 
/(p,f) for a heavy quark in phase space: 


s / (p .<) = 


dj_ 

dt 


( 12 ) 


- collisions 


It is assumed that there is no external force acting on 
the heavy quark and that the phase space distribution / 
does not depend on the position of the quark. The right 
hand side of Eq. represents a collision integral given 

by, 


R(p,t) = 


Of 

dt 


collisions 


= / d 3 k [w(p + k, k)/(p + k) — ru(p, k)/(p)], 


(13) 


as Eq. (|j), earlier. We have also introduced the Bose 
enhancement (Pauli suppression) factors, g(q) = l±g(q) 
for gluons ( quarks) Q. The matrix elements \M\ 2 for 
elastic processes can be obtained from nn- 

We may further simplify the above by assuming soft 
scatterings according to Landau [^| to obtain, 


R 


d 3 k 


i ^ 1 d 2 

k 'dp + 2 kikj d Pi d P j 


(W)- 


(15) 


Now Eq. (ED reduces to well known Fokker-Planck 
form 


Ai{p)f + 3^7 (^'(P)/) 

where the kernels 

A, = j d 3 kw(p, k)ki , 

Bij = i f d 3 kw(p,\i)kikj , 


0J_ = 0_ 

dt dpi 


(16) 


(17) 

(18) 


can be identified with the drag and diffusion coefficients, 
respectively, using the Langevin formalism 1121. In par¬ 
ticular, Aj and are given by, 


1 f d 3 q f d 3 q' r d 3 p' 

1 ~ 2E p J (2n) 3 2E Ci J ( 2 tt ) 3 2 E^ J (2n) 3 2E p , 

x ^-^^^i^f^ip + q-p'-q) 

x 5(q)s(q , )[(p-p')J 

= {{(P-pV) > (19) 

B ij = ^ P)i {p' - P)j)) ■ (20) 


We note that the Eqs. (19) 
and thus we may write [12] 


and (Et«) depend only on p, 


Ai = p.iA(p 2 ) 


Bij — 


A PiPj 

Oij 9 


B 0 (p 2 ) 


PiPj 


Bi(p 2 ) 


( 21 ) 

( 22 ) 

(23) 


where u>(p, k) is the rate of collisions which change the 
momentum of the charmed quark from p to p k. The 
Eq.(|l3l) can be written explicitly Il2) as 


R{pR) = 


1 


d 3 q 


d 3 q' 


d 3 p' 


2E p J (2tt) 3 2E q J (27r) 3 2E q , J (2tt) 3 2E p > 


X — ^ 2 \M\ 2 ( 2 TrfS 4 (p + q-p'-q') 
7c ^ 

x [/(p')3(q')ff(q) - /(p)ff(q)s(q')]» 


(14) 


where 


A = «l»-«p.p'»/p 2 


B 0 = 
Bi = 


i [((p /2 ))-(((p'-p) 2 ))/p 2 ] , 
^[(((P , -P) 2 ))/P 2 -2(( P '-P))+P 2 ((1))]. 


(24) 

(25) 

(26) 


The integrals appearing in the equations given above can 
be further simplified, by solving the kinematics in the 
center of mass frame of the colliding particles, so that we 
can write, 


where p' = p k, q' = q + k, and y c is spin and color 
degeneracy of the charm quarks. In the above, g repre¬ 
sents the distribution for quarks or gluons approximated 
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W_/ (cosx) 

y/(s + ml - m 2 g(q) ) 2 - Asm 2 c 


g( E q) 


/ JL _ /*Z7T 

dcos 0 c . m . Em 1 #«.. 

-1 Jo 

x e^<.'<?(£„,)F(p') . (27) 


where and s = (£ ; p +i? q ) 2 -(p+q) 2 , Eq> = Eq+E p -E p > 
and p' is a function of p, q and 0 C .m.- Note that in ad¬ 
dition to other differences due to the introduction of the 
mass of quarks and gluons (see later) and the quantum 
statistics, this expression differs |2(| by a factor of 2 from 
the corresponding Eq. (3.6) of Ref. Q. The quantum 
statistical correction for gluons employed here introduces 
a divergence as q —> 0 which is absent if we use a Boltz¬ 
mann distribution for them. We have avoided this diver¬ 
gence Q by using the thermal gluon mass; 


m 


2 

a 


— ^ 9 



g 2 T 2 
3 


(28) 


where g is QCD coupling constant (a s = g 2 /An). In order 
to retain a certain degree of self-consistency, we have also 
used the thermal quark mass; 

(»> 


We thus correct these masses for the absence of chemical 
equilibrium. 

We have also evaluated the integrals appearing in 
Eqs. (ED etc. by a direct Monte Carlo procedure, and 
the results were found to agree with the simplified ex¬ 
pressions given above. 


IV. RESULTS 

A. Chemically Equilibrated Plasma 

As a first step we show in Fig. 1(a) the variation of 
the drag coefficient A with the momentum of the charm 
quark p for a chemically equilibrated QGP at a temper¬ 
ature of 500 MeV. First of all, we note that the conse¬ 
quences of introducing the quantum statistics are quite 
small. We must comment on the difference between our 
results and those of Svetitsky before proceeding. 
Firstly, we have used a s = 0.3, instead of 0.6 used by 
Svetitsky. Secondly, we have used Nf, the number of 
flavors, as 2.5 instead of 2. Finally, we have have used 
the Debye screening mass fi = pn = \/Ana s X g T in the 
internal gluon propagator in the t— channel exchange 
diagram, to remain consistent with the description used 
while evaluating the chemical evolution of the QGP. This 
is in addition to the numerical factor of 2 mentioned 


above. We have verified that reverting to the values used 
by Svetitsky, we reproduce the results of Ref. 0 , but 
for the numerical factor of 2. 



FIG. 1. The drag (A) and the diffusion (Bo and B\ — Bo) 
coefficients for charm quarks in a fully equilibrated quark 
gluon plasma at a temperature of 500 MeV. The solid curves 
give the results with the proper quantum statistics while the 
dashed curves give the results with the Boltzmann distribu¬ 
tion. 
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We find that our result for the drag coefficient A is 
a factor about three smaller than obtained by Svetitsky. 
(Note that the results scale with a 2 s .) This will have an 
important consequence as we shall see later. 

The momentum dependence of the diffusion coeffi¬ 
cients B 0 and B\ — B 0 are shown in Figs. 1(b) & 1(c). 
While the general behavior remains similar to the results 
of Svetitsky, our values are smaller, due to the reasons 
given earlier. 



FIG. 2. The energy loss ( dE/dx ) of heavy quarks as a func¬ 
tion of their momenta for T=0.5 GeV and a 3 = 0.3. 

The energy loss dE/dx is related to the drag coefficient 

A; 

TP = -Mp 2 )p • (30) 

dx 

In Fig. 2 we have compared our results for the energy loss 
suffered by charm and bottom quarks with the results of 
Braaten and Thoma jl8| and Bjorken |2lJ (adopted to 
the case of heavy quarks). Our results agree with those 
of Braaten and Thoma within 10%. We note that even 
a fully equilibrated plasma, whose temperature is kept 
fixed at 500 MeV, a charm quark having a momentum 
of 1-5 GeV will travel for almost 8-10 fm before coming 
to rest. The results would be completely different for a 


dE/dx of ~ — 2 GeV/fm assumed in the studies reported 
in Ref. ||IJ,|T5|. We feel however that at exceptionally 
high energies, the additional mechanism of gluon radi¬ 
ation may enhance the energy loss. That should be of 
interest if we wish to study the propagation of a heavy 
quark jet in the QGP. 




FIG. 3. The temperature dependence of the drag (A) and 
the diffusion (Bo and B i — Bo) coefficients for charm quarks 
in a fully equilibrated quark gluon plasma. The results shown 
employ the proper quantum statistics for the gluon and the 
quark distributions. 
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The temperature dependence of the drag and the dif¬ 
fusion coefficients for a fully equilibrated QGP is shown 
in Figs. 3(a-c). We see that the drag coefficient drops 
rapidly with decline in temperature, so that a low tem¬ 
perature QGP hardly offers any resistance to the motion 
of massive quarks, whose final momentum will thus be 
decided by the time it spends in the very hot plasma. 

B. Results for Chemically Equilibrating Plasma at 
RHIC & LHC Energies 

The results given so far are for a fully equilibrated 
plasma. We have already discussed that the plasma likely 
to be created in relativistic heavy ion collisions is far from 
chemical equilibrium. In Fig. 4(a) we give the evolution 
of the temperature and fugacities for the plasma likely 
to be created at RHIC, using the procedure disussed ear¬ 
lier. The appropriate drag and diffusion coefficients at 
any time t operating on the heavy quark are obtained 
by accounting for the current fugacities and the temper¬ 
ature. In order to avoid confusion, we define 

a(p 2 ,r) = A (p 2 , T(t), A 9 (t), A s (t)) , (31) 

and similarly /?o = Bq, Pi = B\. The time dependence 
of the drag and the diffusion coefficients for some typical 
values of the momenta are given in Figs. 4(b-d). We see 
that the drag and the diffusion coefficients are large only 
at the very early times (due to the high temperatures 
then) and drop rapidly as the plasma cools. 

The corresponding results for the plasma created at 
the LHC energies are given in Figs. 5(a-d). We now see 
that the drag and the diffusion coefficients are a factor 
of « 3 larger at any given time for the same momentum 
compared to the case at RHIC. Thus we see that the 
rapid cooling of the plasma and the small values of the 
fugacities ensure that the charm quarks will experience 
strong drag force only very early in their life both at the 
RHIC and the LHC. The force is also unlikely to cause a 
complete stopping and diffusion of the charm quarks. 






t (fm/c) 


FIG. 4. (a)The gluon fugacities (A g ), the light quark fu¬ 
gacities (A g ) and the temperature (T in GeV) as a function of 
proper time (r) for a longitudinally expanding plasma likely to 
be created in Au+Au collisions at RHIC. (b-d) The variation 
of the drag (a) and the diffusion (/3o and (3\ — 0o) coefficients 
with proper time (r) for different charm quark momenta. 
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The final momentum of the charm quark is likely to 
carry the effects of energy loss which in turn is affected 
by the cooling as well as the chemical evolution of the 
plasma. Thus for example, if we had a fully equilibrated 
QGP at the same initial temperatures as used here, then 
the drag at RHIC would be larger by more than a factor 
of three while at LHC energies, it would be large by a 
factor of more than two. Such large values for the drag 
force would ensure that at least the slower charm quarks 
produced at the LHC energies would come to rest soon 
after the creation and diffuse later in a manner discussed 
by Svetitsky. The final momenta of such charm quark (or 
D meson) would then be decided by the hadronization 
temperature. 
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FIG. 5. Same as Fig. 4 for LHC energies. 

In order to get an idea of the energy loss suffered by 
the charm quarks produced at RHIC and LHC, we have 
evaluated their final momentum after they travel through 
the plasma for the duration of the QGP phase. 

Thus we write 

^ = -a(pV)p (32) 

and solve for p as a function of r. The results shown in 
Fig. 6 are very revealing. We find that by the time the 
QGP phase is over, the charm quark which was produced 
at the beginning of the collision would have lost up to 
40% of its initial momentum at LHC energies, while at 
the RHIC, the energy loss may not exceed 10-12%. 











FIG. 6. The evolution of the charm quark momentum as a 
function of time in a partially equilibrated plasma at RHIC 
and LHC. The charm quarks are assumed to be produced with 
momenta of 1, 2, and 3 GeV at the time to. 


C. Distribution of Charm Quarks at RHIC & LHC 


This has very interesting consequences. At RHIC en¬ 
ergies, one may thus treat charm quark, almost as pene¬ 
trating probes, which will provide information about the 
momentum distribution of the initially produced charm 
quark pairs. One can estimate the momentum distri¬ 
bution of initially produced charm pairs in a central 
Au+Au collisions by scaling the heavy-quark pair pro¬ 
duction cross-section in pp collisions as p8fl 


dN 

d 2 pTdy 


—TabQ) = 0 ) 

7T 


da 

dp^dy 


(33) 


where Tab is the nuclear thickness factor. For central 
Au+Au collisions it’s value is 293.2 fm -2 ^9|. One may 
further approximate the thermal production of the charm 
quark pairs from Eq.(25) of Ref. |3(|. The result of this 
study is shown in Fig. 7. One may obtain the “energy-loss 
corrected” distribution by decreasing the momenta by 
about 10%. This also implies that the dileptons from the 
annihilation of quarks would continue to remain buried 
under the leptons originating from the open charm de¬ 
cay, as originally inferred by authors of Ref. jll] . We may 
also add that as the charm quarks do not stop/diffuse at 
RHIC, they will not be affected by the transverse veloc¬ 
ities which may develop during the QGP phase. 

The situation is more complex at the LHC energies. 
Now the charm quarks can not be treated as penetrating 
probes. The results given in Fig. 8 correspond to situa¬ 
tion when we treat them as penetrating probes, as in the 
early studies in this field. As most of the heavy quarks 
are produced very early in the collision, we can still get 
an idea of the final momentum distribution by decreasing 
the momenta in this figure by about 40%. This however 
may not be quite appropriate, as the charm quarks which 


have lost a large fraction of their initial momentum due to 
the drag force will start getting affected by the transverse 
flow of the QGP which can grow to large values at the 
LHC energies |j|. It is not clear that this situation can 
be accurately handled within a hydrodynamic descrip¬ 
tion of the expansion. A more appropriate description 
could be the parton cascade model which includes colli¬ 
sions and even radiations from the charm quark. This 
work is under progress. 



Pt (GeV) 


FIG. 7. The momentum distribution of charm quarks for 
RHIC energies without energy loss. 



Pt (GeV) 


FIG. 8. The momentum distribution of charm quarks for 
LHC energies without energy loss. 

All the same it is not difficult to imagine that the de¬ 
crease in the slope of the pt distribution of the charm 
quarks is offset to a large extent by the flow velocity likely 
to develop at the LHC energies, and thus the dileptons 
from the plasma will perhaps remain buried under the 
open charm decay, at these energies as well, unless we 
look at too large masses. 
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V. SUMMARY 

In brief, we have obtained the drag and diffusion co¬ 
efficients for charm quarks in a chemically equilibrating 
quark gluon plasma which may be produced at RHIC 
and LHC energies. Using a set of reasonable initial con¬ 
ditions, we find that a charm quark produced at the early 
stage of the collision may loose up to 10% of its momen¬ 
tum during the life-time of the QGP phase at the RHIC 
energies. This suggests that at the RHIC energies, dilep¬ 
tons originating from annihilation of quarks may remain 
buried under the background from open charm decay. 
The situation at LHC is more complex, as the charm 
quarks may loose up to 40% of the initial momentum. 
This could however be somewhat offset by the transverse 
flow of the QGP fluid, which could be large there. Fur¬ 
ther work, preferably within a parton cascade model is 
needed to settle some of these questions in a more definite 
manner, though we feel that even at the LHC energies, 
the dilepton signal may remain buried under the back¬ 
ground from open charm decay. We may add that some 
recent works in this direction have used too large values 
for the energy loss, and arrived at results differing from 
our observations. 

The small energy loss seen for heavy quarks in an equi¬ 
librating plasma could lead to a difference in the “quench¬ 
ing” of a heavy quark jet as compared to that of a gluonic 
jet. (Radiative loss of gluonic jets will presumably not 
be reduced, considerably, in an unequilibrated plasma.) 
This can be of great interest |jlj . 
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